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The Color Octet Effect from e+e" ^ J/ip + X + at B Factory 

Jian-Xiong Wan^ 
Institute of High Energy Physics, Academia Sinica, Beijing 100039, CHIN.^ 
(Dated: February 1, 2008) 

We study the initial state radiation process e+e~ — > J/-!/' + X + 7 for J/^p production at B factory, 
and find the cross section is 6f % larger than it's Born one for color octet part and is about half as 
it's Born one for color singlet part. Furthermore, the color singlet and color octet signal are very 
clearly separated in it's E-y spectra due to kinematics difference. We suggest to measure this E-y 
spectra at B factory to determine the color octet effect. 
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To study J production on various experiments is 
a very interesting topic since its discovery in 1974. 
It is a good place to probe of both perturbative and 
nonperturbative aspects of QCD dynamics. To de- 
scribe the huge discrepancy of the high-p^ J 1^ pro- 
duction between the theoretical calculation based on 
color singlet mechanism^ and the experimental mea- 
surement by the CDF collaboration at the TevatronQ, 
color-octet mechanismQ was proposed based on non- 
relativistic QCD(NRQCD)||. The factorization formal- 
ism of NRQCD provides a theoretical framework to the 
treatment of heavy-quarkonium production. The color 
singlet mechanism is straightforward from the pertur- 
bative QCD, but the color-octet mechanism depends on 
nonperturbative universal NRQCD matrix elements. So 
various efforts have been made to confirm this mecha- 
nism, or to fix the magnitudes of the universal NRQCD 
matrix elements. Although it seems to show qualitative 
agreements with experimental data, there are certain dif- 
ficulties in the quantitative estimate of the color-octet 
contribution in J/tj) and V'' photoproduction at HERA 
[EIEI^IEIS; J/^ii'') polarization in large transverse mo- 
mentum production at the Fermilab Tevatron [lOlllllIT^ . 
and more recently in B-factories. 

There are at least two kinds of uncertainty in the the- 
oretical treatment of J/ip production. One is from the 
QCD correction. It is difficulty to estimate the QCD cor- 
rection effect up to a few times in various leading order 
calculation for Jf-ip production. The uncertainty to fix 
the NRQCD matrix elements is quit large for the case 
in which color singlet effect is compatible with the color 
octet one such as in the measurement of the J/^ produc- 
tion at B factory in BaBar and Belle experiments Il4l | . 
There are a few examples shown that the NLO correction 
are quit large. It was found that the current experimen- 
tal results on inelastic J/ip phtonproductionp^a. Ilq are 
adequately described by the color singlet channel alone 
once higher-order QCD corrections are included^- Al- 
though ref[l3 found that the new DELPHI ^ data ev- 
idently favor the NRQCD formalism for J/4> production 






FIG. 1: The typical Feynman diagrams for the e"'"e~ —f 
J/tpX, there are 6 diagrams for the e^e~ —^ J/i^gg, 4 di- 
agrams for e'^e' — > J/ipcc and 2 diagrams for each e^e~ 
^pP|]g, J = 0,1, 2 



7 -I- 7 ^ J/ip + X, but rather the color-singlet model. 
It was also found by ref 19] that the QCD higher order 
process 7 -f 7 — s- J/ipcc gives the same order and even 
bigger contribution at large compare to leading order 
color singlet processes. refj2flj found that higher order 
process give large contributions. 

The another uncertainty is from the hadronization 
of the color octet J/-)/)'*' to the color singlet J/ip and 
hadrons. The recent measurement of the J/tjj produc- 



tion at B factory in BaBar and Belle experiments fisl 
show that the theoretical predication pH for Pj/^f, 
spectra did not agree with experimental results, and 
gives the total cross section a{e~^e~ J/tp + X) = 
1.47ib0.10ib0.13p6. It seems compatible with theoretical 
prediction 23. 2112511 cover the range 0.8— 1.7 pb , in which 



about O.Spb from color singlet processes e+e 
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^ J/tpcc, and 0.5 — 0.8pb from color octet 
processes e+e~ J/ip^^^g. The color octet processes 
only contribute to the endpoint of p}/^ spectra due to 
the kinematics of the two body final state. But the ex- 
periments did not observe the this signal. To explain 
this discrepancy, it is nature to think that J /ip^^'> have 
to hadronize into color singlet J /'ip and will lose it's en- 
ergy such as the case when a quark jet hadronize into 
hadrons. Therefore there should be a hadronization pos- 
sibility function Fix,...) for J/%p production with mo- 
mentum pj/^ — X Pj/^(s) and Jq dx F{x,...) — 1. The 
universal NRQCD matrix elements treatment is just a 
first step approximation oi F{x,...) = 5{1 — x). The 
refp^ tried on it for the J/^ production at B factory 
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FIG. 2; The typical Feynman diagrams for the e e~ —* 
J/'!^X7, there are 12 diagrams for the e^e~ — > J/ipgg^, 
8 diagrams for e^e~ J/iJjcc'^ and 4 diagrams for each 
e+e- ^pPj]g7, J = 0,1, 2 

and broadened the p* spectra from color octet J /ip- 
Due to the above two kinds of uncertainty, the Pj/^ 

spectra from color octet J/i/''-*^ is unknown, and mean- 
while the higher order QCD effect for the color singlet 
J production is also unknown. So to confirm, to ex- 
tract J/t/j^^'^ information from the measurement in B fac- 
tory experiment become unavailable. However, we find a 
good way to determine the color singlet and color octet 
effect, i.e. to measure the spectra in the initial state 
radiation process e^e^ — > J/ip + X + which could 
avoid the uncertainties from J/tp^^'^ hadronization and 
QCD correction. 

This QED higher order process is thought to be small 
at first glance. But there are already few examples in 
which QED higher order process could be very lager in 
certain case. It was found that there are large contri- 
butions to J/V' production form QED processes such as 
e'^e" — > J/'07 and e"*"e~ — s- Jl'ipe^e~ 27] due to en- 
hancement from the t-channel peak. In the try to ex- 
plain the double J/ij} production measured by Belle |[28l|. 
refp^ found that the high order QED process e+e" 
J/ip + J/ip gives large contribution. Our numerical re- 
sults show that e+e~ — > J/ipX + 7 gives larger cross 
section than e^e~ J/4> + X . 

For the leading process e+e^ J/^ + X, the typ- 
ical Feynman diagrams are shown in Fig^ The color 
singlet processes e~^e~ J/ipgg, J/ijjcc. and color octet 
processes e+e~ — > ^ pPT j q, J = 0,1,2; V[^5'o]5 were 
calculated in ref.^J, iS Hi El- For the e+e" 
J /ip -\- X + 7, there are e"'"e~ — *■ J/ipggi, J/ipcc + 7, 
i^[^Sl]g-i, ^[^Pj]gi, J = 0,1,2, and the typical Feyn- 
man diagrams are shown in Fig|21 

To calculate this initial state QED radiation correc- 
tion, we have to include two parts. One is the virtual 
correction plus soft photon emitted process, the another 
is the hard photon emitted process. A universal formula 
is easily obtained for the first part, i.e. virtual correction 
plus soft photon emitted one. It isfsoj: 

a^"^ = ao{ — [(In — - l)(ln^ + j) + IT " A^^' 
TT y/s 4 5 4 

Where the ctq is the Born cross section, the Se is maxi- 
mum energy of the emitted soft photon, s is the center 



TABLE I: The numerical results for the color singlet and 
color-octet J/t/j production. The Pj/^ = — ^'^J/i> 

the 3-momentum of J/ip. Se = 0.1 is used. The cut 
Pj/^ > 2.0GeV is used by the experiment to suppress the J/ip 
from B decay. The result for color octet with cut condition 
p*j/^ > 2.0Gey is not presented since it does not make sense 

in the case that the detail of J/?/)'*' hadronization procedure 
is unknown 





Born 


VirtuaH-Soft 


hard photon 


p*j^^ > 2.0GeV 


e+e -> J/ipgg{"f) 


0.158 


-0.044 


0.072 


e^e~ — » J/ipcclj) 


0.102 


-0.028 


0.022 


In total 


0.270 


-0.072 


0.094 


pj/^ > OGeV 


e+e ^ J/^iiggin) 


0.203 


-0.056 


0.099 


e+e~ — » J/ij)cc{')) 


0.119 


-0.033 


0.029 


In total 


0.322 


-0.089 


0.128 


e+e- ^ ^rP5]fl(7) 


0.143 


-0.039 


0.351 


e+e- ^ 4,['Sl]g{^) 


0.133 


-0.037 


0.094 


In total 


0.276 


-0.076 


0.445 



mass energy of e^e~ collider. In the above formula, the 
infrared divergence are canceled between the soft photon 
emitted process and virtual correction. 

For both the Born processes and the hard photon 
emitted processes, We obtain the formula and Fortran 
source of the cross section, including the kinematics, by 
using our automatic Feynman diagram calculation pro- 
gram FDC 31]. The FDC was used in the calculation of 
e+e- j'/ip-i and e+e" Jj^e+e' [13]. In refpl] it 
was used in the calculation of 7 -I- 7 —> J/ V'cc, and many 
other processes to repeat the calculation by ref0]. In 
the Fortran source generated by FDC, it is very easily to 
check the gauge invariance by replacing the polarization 
vector of gluon or 7 to it's momentum. All the possible 
gauge invariance are checked to be satisfied numerically 
in our calculation. Our analytic formula generated by 
FDC for e+e- ^ V'P^'jlff, J = 0,1,2; ^'^Sl\g are 
checked with ref^2lj. 

The numerical results are listed in Tableland shown 
in Fig|3| and ^ In the numerical calculation, we use 
nic — 1.5, rriji^ = 3.0, a — 1/137 and the lowest-order 

formula for ai"^"^'' (/i). J jili matrix elements are chosen 
as < J/V' >= 1.4. The universal NRQCD matrix ele- 
ments for color octet J j^ip are chosen as that in ref0],i.e, 

< Jl^^^^{^Sx) >= 0.0039, < J/?/'(^H^5o) >= 0.015, 

< tp^^^^Pj) >= (2J-I- l)0.0043mc2. 

The renormalization scale /i dependent in the calcula- 
tion are as = ai'^^(10.6) — 0.188 for Born processes and 

as — a'f\^y (pi +P2 — P-yY) for hard photon processes. 
Where the pi , p2 and p^ are the momentums for e+ , e~ 
and the emitted 7 respectively. The hard photon pro- 
cesses are defined as > Se , i.e. the 6e is the minimum 
energy of the emitted photon. We have used different 
value of Se = 0.5, 0.4, 0.3, 0.2, 0.1, O.OlGeF to verify that 
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FIG. 3: The differential cross section distribution vs the en- 
ergy of emitted hard photon 
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FIG. 4: The differential cross section distribution vs the 
cos{9) of the emitted hard photon and the beam 



the Eq^ which included only the leading term from 
the soft photon, is a good approximation when 5e is small 
enough. The numerical results show that the result of vir- 
tual -I- soft + hard is almost independent of <5e. There is 
infrared divergence for the process e^e^ '4'^Pj]9{l) 
when the energy of gluon approach to zero, but it can 
not happen since the gluon must hadronize into at least 
one pion. To apply the same cut condition used by the 
experiment ^4j| to suppress the background, the condi- 
tion \/ (pi + P2 — P-yY > "^J/ip + 3m7r is applied to all 
the radiation processes. 

There are other processes which may give compatible 
contributions from the perturbative order analysis, but 
the contributions form them are found too small so that 
they do not appeared in the above list. Such as e+e~ 
V'p5'f]57 is 1.2 X 10" V- e+e- J/ipgge+e- with 48 
Feynman diagrams are calculated and the result is 6.4 x 

In Fig01 we find that most of the emitted photon go 
out with small angle to beam direction in both the color 
singlet and color octet processes. This is the t-channcl 
enhancement effect which make the contribution of the 
QED correction large. In FiglSJ it very clearly shown 
that the color singlet processes mainly emit the photon 
with low energy [E^ < l.OGeV), meanwhile, the color 
octet processes mainly emit the photon with high energy 
{Ej > A.OGeV). The reason for this distribution is com- 
pletely form the kinematics structure as shown in FigI21 
When the center mass energy go down from 10.6Gev to 
5GeV, the cross section will increase about 4 times for 
e+e" J/ipgg, about 51 times for e+e" ''P[^Pj]9 and 
about 17 times for e+e" — *■ i(j[^S'~']g, and the channel 
e'^e~ —> J/ipcc will be closed. The above property was 
shown in the Fig. 5 in ref(23j and the Fig. 3 in ref[25l|. 



So that the remaining subprocesses will show this prop- 
erty after the photon emitting from the initial electron 
or positron. The peak in the small E^ is well known 
from the soft photon effect. Based on the above analy- 
sis, it is very clear that the spectra in FigOl comes from 
the kinematics structure of the hard photon emitted pro- 
cesses. It is hardly be changed even if there is large QCD 
correction effects since the QCD correction appears only 
in the subprocesses and will not change the kinematics 
structure. So the ratio of the cross section for hard pho- 
ton process to that for it's Born processes will be hardly 
changed by the QCD correction. In another side, the E^ 
spectra can not be changed by the procedure of J/tp^^'^ 
and gluon hadronization into J/ip and hadrons since the 
hadronization procedure has no relation with it. 

In the measurement of J/iJj production at B factory in 
BaBar and Belle experiments [l3 . Il3 | . it was tried to sepa- 
rate the color singlet and color octet effects by measuring 
the spectra oi p*j^^. The spectra show that the theoret- 
ical predication for color octet does not appear. It let 
us have to think about the detail of how the J/V'^*-' go 
through hadronization procedure. The ref[26| addressed 
on the problem, but the results is not predictive since 
the introduced shape function for J/^'^^ is unknown. To 
avoid the uncentainties from J/f/j^^-* hadronization and 
QCD correction, we find that to measure the E^ spec- 
tra of the hard photon emitted J/ip production process 
e^e~ — !■ J/V' -l- X -|- 7 in B factory experiment is a very 
clear and good way. This spectra is not changed by QCD 
correction, by the detail of J/ip'-^'> hadronization. and 
will very clearly separate the color singlet and color octet 
signal. Another advantage is that the color octet cross 
section in the hard photon process is even 61% large than 
it's Born one. 
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Furthermore, It has been argued in the above that the 
ratio of the cross section of the hard photon process to it's 
Born one will be hardly changed by QCD correction and 
the detail of J/ip^^'> hadronization. As approximation, 
we have 

a{e+e- ^ J/V' + X) ^ 0.233a; + 0.2y, 
cr(e+e- ^ J/^ + X + j) ^ 0.128X + 0.445y, 
cr(e+e^ — s- J/^jj + X, color singlet) = 0.233a;, 
cr(e^e^ — > J/ip + X, color octet) — 0.2?/. 

where x and y represent the effects of QCD correction 
and uncertainty from the NRQCD matrix elements, and 
all the numerical values from the TableQ] The color octet 
and color singlet contributions can be extracted by solv- 
ing EqElwhen both the cross section cr(e+e^ J/ip+X) 
and ^(e+e" — > J/ip + X + 7) be measured in the B fac- 
tory experiment. In this way, the color octet cross sec- 



tion could not apply cut condition P*j/^ > 2.0GeV since 
the P*j/^ spectra depends on J/ip^^^ hadronization. So 
how to suppress the background from B decay becomes 
a problem. However the measurement of the E-y spectra 
will not suffer from this problem and the E-y spectra fit- 
ting will fix the contributions from color singlet and color 
octet mechanism. 



Acknowledgments 

This work was supported in part by the National 
Natural Science foundation of China under Grant 
Nos. 90103013 and by the Chinese Academy of Sciences 
under Project No. KJCX2-SW-N02. 



[1] C.H. Chang, Nucl. Phys. B172 (1980) 425; 

J.H. Kiihii, J. Kaplan, and E.G.O. Safiani, Nucl. Phys. 
B157 (1979) 125; 

B. Guberina, J.H. Kxihn, R.D. Peccei, and R. Riickl, 
Nucl. Phys. B174 (1980) 317; 

E.L. Berger and D. Jones, Phys. Rev. D23 (1981) 1521; 

R. Baier and R. Riickl, Z. Phys. C19 (1983) 251. 
[2] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 69 

(1992) 3704. Phys. Rev. Lett. 79 (1997) 572. Phys. Rev. 

Lett. 79 (1997) 578. 
[3] E. Braaten and S. Fleming, Phys. Rev. Lett. 74 (1995) 

3327 

[4] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. 

D 51 (1995) 1125 [Erratum-ibid. D 55 (1997) 5853] 
[5] M. Cacciari and M. Kramer, Phys. Rev. Lett. 76 (1996) 

4128. 

[6] J. Amundson, S. Fleming, and I. Maksymyk, Phys. Rev. 

D56 (1997) 5844. 
[7] P. Ko, J. Lee, and H.S. Song, Phys. Rev. D54 (1996) 

4312; [60 (1999) 119902(E)]. 
[8] B.A. Kniehl and G. Kramer, Phys. Lett. B413 (1997) 

416. 

[9] M. Kramer, Nucl. Phys. B459 (1996) 3. M. Kramer, J. 
Zunft, J. Steegborn and P. M. Zerwas, Phys. Lett. B348 
(1995). 

[10] M. Beneke and LZ. Rothstein, Phys. Lett. B372 (1996) 
157 [Erratum-ibid. B389 (1996) 769]; M. Beneke and M. 
Kramer, Phys. Rev. D55 (1997) 5269. 

[11] E. Braaten, B.A. Kniehl, and J. Lee, Phys. Rev. D62 
(2000) 094005; B. A. Kniehl and J. Lee, Phys. Rev. D62 
(2000) 114027. 

[12] Adam K. Leibovich, Phys. Rev. D56 (1997) 4412. 

[13] B. Aubert et al. [BABAR Collaboration], Phys. Rev. 
Lett. 87 (2001) 162002 

[14] K. Abe et al. [BELLE Collaboration], Phys. Rev. Lett. 
88 (2002) 052001 

[15] HI Collaboration, Contributed Paper 157aj, Interna- 
tional Europhysics Conference on High Energy Physics 
(EPS99), Tampere, Finland, 1999. 

[16] ZEUS Collaboration, Contributed Paper 851, Interna- 



tional Conference on High Energy Physics (ICHEP2000), 
Osaka, Japan, 2000. 

[17] M. Klasen, B.A. Kniehl, L.N. Mihaila, and M. Stein- 
hauser, Phys. Rev. Lett. 89 (2002) 032001. 

[18] S. Todorova-Nova, in Proceedings of the 31st Interna- 
tional Symposium on Multiparticle Dynamics, Datong, 
China, 2001 [ArXiv: hep-ph/0112050 ; M. Chapkin, talk 
at 9th International Conference on Hadron Spectroscopy, 
Protvino, Russia, 2001, AIP Conf. Proc. 619 (2002) 803, 
hep-ph/0307049;to appear in PLB). 

[19] C. F. Qiao and J. X. Wang, to be appear in Phys. Rev. 
D. [ArXiv: hep-ph/0308244 

[20] C. F. Qiao, J. Phys. G: Nucl. Part. Phys. 29 (2003) 1075 

[21] E. Braaten and Y. Q. chen, Phys. Rev. Lett. 76 (1996) 
730 

[22] Peter Cho, Adam K. Leibovich, Phys. Rev. D 53 (1996) 

6203 Phys. Rev. D 54 (1996) 6690 
[23] F. Yuan, C.-F. Qiao, and K.T. Chao, Phys. Rev. D 56, 

(1997) 321. K.T. Chao and L. K. Hao, Nucl. Phys. B115 

(2003) 162. 

[24] S. Baek, P. Ko, J. Lee and H. S. Song, J. Korean Phys. 

Soc. 33 (1998) 97 
[25] G. A. Schuler. Eur. Phys. J. C8,273(1999) 
[26] S. Fleming, A. K. Leibovich and T. Mehen. 

arXiv:hep-ph/0306139 . 
[27] C.H. Chang, C.F. Qiao, and J.X. Wang, Phys. Rev. D56 

(1997) 1363; 57 (1998) 4035. 
[28] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 89, 

(2002) 142001. 

[29] G. T. Bodwin, J. Lee and E. Braaten Phys. Rev. Lett. 

90 (2003) 162001 
[30] A. B. Arbuzov, V. A. Astakhov, A. V. Fedorov, G. V. Fe- 

dotovich and E. A. Kuraev JHEP 9710:006,1997, [ArXiv: 

hep-ph/9703456 
[31] J. X. Wang, FDC project was started at 1992 and 

many new parts added since then. "Proceeding of AI93, 

Oberamergau, Germany. New computing techniques in 

physics research HI" 517-522," Appi 2001, Accelerator 

and particle physics" 108-121. 



